Macroscopic Greenberger-Horne-Zeilinger and W States in Flux Qubits 
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We investigate two types of genuine three-qubit entanglement, known as tlie Greenberger-Horne- 
Zeilinger(GHZ) and W states, in a macroscopic quantum system. Superconducting flux qubits are 
considered theoretically in order to generate such states. A phase coupling is proposed to offer 
enough strength of interactions between qubits. While an excited state can be the W state, the 
GHZ state is formed at the ground state of the three flux qubits. The GHZ and W states are shown 
to be robust against external flux fluctuations for feasible experimental realizations. 

PACS numbers: 74.50.+r, 85.25.Cp, 03.67.-a 



Introduction.— FiTitaiiglemeTit plays a crucial role in 
quantum information science. Controllable quantum sys- 
tems such as photons, atoms, and ions have provided 
the opportunities to generate the entanglements. Re- 
cent experiments on two qubits have shown the exis- 
tence of entanglement in different types of microscopic 
systems. Further, multipartite entanglements such as 
the Greenberger-Horne-Zeilinger(GHZ) [1] and W 0, Q 
states have been demonstrated in recent experiments of 
atoms 01, photons and trapped ions [liQ- But, in solid- 
state qubits it has not yet been achieved. 

As a macroscopic quantum system, superconducting 
qubit systems have been investigated intensively in ex- 
periments because their system parameters can be con- 
trolled to manipulate quantum states coherently. In- 
deed, the entanglements between two charge 0, phase 
[liQ, and flux qubits 1^, 11 1 have been reported. While 
the timely evolving states in the experiments of charge 
qubits exhibit a partial entanglement, the excited level 
(eigenstate) of capacitively coupled two phase qubits Q 
shows higher fidelity for the entanglement. The experi- 
ments in Ref. [lo| show a possibility that two flux qubits 
can be entangled by a macroscopic quantum tunneling 
between two-qubit states, flipping both qubits. Actu- 
ally, the higher fidelity in the capacitively coupled two 
phase qubits is caused by the two-qubit tunneling pro- 
cesses [9|. In a very recent study, the two-qubit tunnel- 
ing process was theoretically shown to play an important 
role in generating the Bell states, maximally entangled, 
in the ground and excited states |13 |. 

For multipartite entanglements in superconducting 
qubit systems, there have been few studies. To produce 
the GHZ state in three charge qubits, only a way of doing 
a local qubit operation via time evolutions was suggested 
[l^. As one of possible directions to produce such mul- 
tipartite entanglements, then, it is natural to ask how 
to create the W state as well as the GHZ state in the 
eigenstates of superconducting three-qubit systems. Here 
we consider three flux qubits. Normally, the interaction 



strength between inductively coupled flux qubits [l^ is 
not so strong that the controllable range of interaction is 
not sufficiently wide. To control a wide range of interac- 
tion strengths in the qubits, we use the phase-coupling 
scheme [11 H El, H, [l?! for three qubit (see Fig. [T]^a)) 
which enables to generate the GHZ and W states and to 
keep them robust against external flux fluctuations for 
feasible experimental realizations. 

Model.— We start with the model shown in Fig. ^l[a)■ 
The Hamiltonian is written hy H = 

Ucs{4>), where Pi = —ihd/dcpi and My = 



2 I ij 



P. 



($o/27r)2C,% 

with the capacitance of the Josephson junctions Ci. 
The dynamics of the flux qubits [20| are described by 
the phase variables = {cjjqi, (j)'^) with q = a,b,c and 
i = 1,2,3, where (p^s are the phase differences across the 
Josephson junctions. If we neglect the small inductive 
energy, the effective potential is written in terms of the 
Josephson junction energies, Ucs{(t>) = SgELi 
cos (j)qi) + E'j{l — cos(j)'q)]. The periodic boundary con- 
ditions involved in the qubit loops and the connecting 
loops can be written as 

(1) 
(2) 
(3) 



i>ql + 4>q2 + 093 = 27r(71g + fq), 

(0ai - 0ci) - 0a + 0c = 27rr, 

(0fcl - 0cl) - 0fc + 0c = 27rs, 



where q = a,h,ci& qubit index and r, s, riq integers. Here 
fq = <f>g/$o with external flux <I>g and the superconduct- 
ing unit flux quantum <f>o = h/2e. Two independent 
conditions in Eqs. ([2]) and ([3]) are the boundary con- 
ditions for connecting loops. For simplicity we consider 
E,j2 — Ej3 — Ej and C2 = C3 = C, so we can set 
092 = 093 and Eq. ^ becomes 0,1 = 2t: {uq + fq) ~2(j)q3. 
The results for Ej2 ^ Ej^ are qualitatively the same. 

At the coresonance point {fa, fb, fc) = (1/2, 1/2, 1/2), 
the effective potential is given by 

f/eff(0) 



(-Bji COS 20^3 
I- 3Eji+6Ej + 3E'j 



2Ej cos 0,3 - E'j cos (j)q) 
(4) 
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Here, we introduce a rotated coordinates (p = 
{4>aT4>i3^4>i) in Fig- Illlb). The Euler rota- 
tions provide new coordinates such as (p^ = 

7^2(x,_o,o)7^l(o,o,0)¥'^ = 7^(x,e)<P^ with x = 

-tan"iV2, 9 = -7r/4 and ip = ((/>a3, '/'bs, ^cs), which 
can be written explicitly as 



1 



V3 -V3 
^ , 1 1 -2 
^6 \ 72 ^/2 V2 







1 






Uc3/ 



(5) 



In the same way, a new coordinates for ip' = {(j)'^, cji'^^, cj)'^) 
is given as ^'^ = 7^(x, e)ip'T with ^' = ((/)^, 0^, Us- 
ing the boundary conditions of Eqs. (HI)-® the Hamil- 
tonian is written in the transformed coordinates, (p = 

(0a, 0/3, '^'7)' as = E^=a,/3,7 + 2M; + t^ofr(^), 

where M« = = 4Ci 2C -I- 2C", = 4Ci -h 2C, and 
C . Note that the value of 6' is determined at the 



potential minimum, dUcff{4')/d4>j = 0. 

The eight corners of the hexahedron in Fig. [Ijb) 
correspond to the three-qubit states. Here the | |) 
(I T)) is defined as diamagnetic (paramagnetic) cur- 
rent state which corresponds to positive (negative) 
value of (j)qi in the boundary condition of Eq. p]). 
These states can be represented more clearly in the 
rotated coordinates, (0ct,0^,0^), because the effective 
potential Ucs{<P) has three-fold rotational symmetry 
about the 0^-axis, which can be shown as in the 
follows. Using the transformation of Eq. ^ one of 
the terms in Eq. (jH) is written as ^q^a b c'^'^^'t'qs ~ 

^ cos(200/V6)] - 
- sin(200/V6)]. 
Here, if we rotate the potential by 2tt/3 about the 0^ axis 
as^c. = -(l/2)0,-(V3/2)00,V/3 = (\/3/2)0„- (1/2)0/3 
and "07 — 07, we can easily check the invariance of the 
effective potential, Uca{4>)- 

In order to study the GHZ state, I^'ghz) = 
(ITTT) + liii)) /V2, we draw the yellow (dark gray) 
square introducing the auxiliary coordinates defined by 

0p = (0a3 + 0h3)/A/2 and 0m = (0a3 - 063)7 ^2, whilc 

for W state, l^-w) = (ITii) + UTi) + liit)) /V3, we con- 
sider the blue (light gray) triangle. Figures [Hc)-(e) show 
the effective potential Ucs{4>) in Eq. ffl. When the 

Fig. He) 
is dominant 
As E'j increases, it 



cos(0^/ V3) [2 cos (00 / %/6) cos(0Q,/ V2) 
sin(0-y/\/3) [2 sin(00/V6) cos(0q/-\/2) 



mu. 



three qubits are decoupled for Ej = 
shows that the single-qubit tunneling, 
over the three-qubit tunneling, 
is shown in Fig. [Ijd) that the three-qubit tunneling be- 
comes dominant. Then the GHZ state is expected to be 
formed at the ground state. The dotted line in Fig. [T]^d) 
coincides with 0-y axis in Fig. [Ijb). Along the 0-y axis, 
the double- well potential is given by Uee{0, 0, 07, VStt) — 

3E,n (l + cos ^) + 6Ej (l - cos ^) , where the bar- 
rier hight is proportional to Eji. The WKB approxima- 
tion allows us to calculate the three- qubi t tunneling, t^, 



through this double- well potential [21|, . Other tunnef 
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FIG. 1: (Color online.) (a) A three flux qubit system. The 
black squares are the Josephson junctions. The Josephson 
coupling energy of the Josephson junctions in the qubit and 
connecting loop are Eji and E'j, respectively, /q's are the 
external fluxes and 0's are the phase differences across the 
junctions, (b) The eight states of three qubits are represented 
in (0a3, 0b3, 0c3)-space at the coresonance point, ((pa, 0/3, 
are the rotated coordinates and 0(0,0,0) is the origin of both 
coordinates. The blue (light gray) triangle intersects verti- 
cally the 0-y axis at point P. For Eji/Ej = 0.7, the effective 
potentials UeS in (0p, 0c3)-plane ( yellow (dark gray) square 
in (b)) are plotted for (c) Ej^O and (d) E'j/Ej = 0.5. 
The dotted line in (d) coincides with 0^ axis in (b). (e) The 
effective potential ?7e£F in {4>a, 0/3)-plane (blue (light gray) tri- 
angle in (b)) for E'j/Ej = 0.05 and Ej-^/Ej = 0.75. Here 
and after, the superscript a{b) in f"'*"' denotes the tunnelling 
processes including (excluding) the states, | HT) or | Hi), 
and i = 1,2,3 the single-, two-, and three-qubit tunnelling 
processes, respectively. 



ings such as single-qubit tunnelings, t° and t\, and two- 
qubit tunnelings, and can also be calculated. The 
tight-binding approximation based on the eight states 
of three qubits gives the effective Hamiltonian, H = 
J2„E'^W)H - Yji,yi'^'^'W)W\^ where t^,,, = tl^''^ and 

\v) = \SaSbSc) with Sq G {T,i}- 

Q-measure.— The global entanglement for tripartite 
systems can be quantified by the Q-measure fisj. For a 
normalized arbitrary three-qubit state, j^*) = ci + 

C2 UiT)+c3 \in)+ci UTT)+C5 \UI)+C6 ITiT)+c7 ITTi> + 

C8 ITTT), the Q-factor is given by 



q(i*)) = |E^.'-(|*))' 

7 = 1 



(6) 



where Di{\^)) = jcice - C2C5P -I- |ciC7 - C3C5P + \cics - 

C4C5P -I- |C2C7 - CsCeP -I- |C2C8 - C4C6p -|- |C3C8 - C4C7I' 
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FIG. 2: (Color online.) The Q-factors of the ground state 
in the three qubit system for (a) E'j/Ej = 0.02 and for (b) 
E'j/Ej = 0.05. Here, = ^3/2 and E.,i/E., = 0.7. (c) 
Cut view of Q-factors in (a) and (b) for fa = 0. (d) For 
fa = fp = and E'j/Ej = 0.6, Q-factors are plotted as a 
function of f-/ for several Eji. 



and 1)2(1 5')) and Dad^t')) are obtained by exchanging 
the indices as 3 5, 4 <-> 6 for D2 and 2 <-> 5, 4 ^ 7 for 
D3. For the GHZ state, Qd^cHz)) = 1 and for the W 
state Q(|*w)) = 8/9. 

GHZ state.— We plot the Q-factors for the ground 
state in Figs. EJa) and (b) as a function of the ro- 
tated fluxes fa = {fa - fb)/V2Jl3 = {fa + fb - 

2/c)/a/6, A = {fa + fb + fc)/VS- Note that the cores- 
onance point {fa, ft, fc) — (1/2, 1/2, 1/2) is transformed 
to {fcjpj^) = (0,0, a/3/2). For E'j/Ej = 0.02 in Fig. 
[2] (a), Q < 1. But, as E'j increases, the GHZ state ap- 
pears around the coresonance point in Fig. [2] (b). Fig- 
ure m (c) is the cut view of Q-factor for various coupling 
strength. It is found that for the GHZ state, Ej should 
be larger than 0.03i?j. It turns out that the couphng 
strength from the inductive coupling scheme corresponds 
to Ej « 0.005£^j [31. This shows that the inductive cou- 
pling scheme cannot provide a sufficient coupling for the 
GHZ state. 

To be observed experimentally, the GHZ state should 
be robust against fluctuations of external flux. Figure 
HKb) shows that the GHZ state can be obtained for a 
broad range of fa and /^. Thus, let us examine the be- 
havior of Q-factor as a function of (Fig. Hfd)). If the 
peak width is too narrow compared with the fluctuations 
of external flux, the GHZ state cannot be observed ex- 
perimentally. Actually, it is found that the three-qubit 
tunneling t^ plays an important role for wide peak width. 
If other tunneling processes except t^ are negligible, small 
flux fluctuations can influence the Q-factor given ap- 
proximately by Qd*)) « {t^)^/{{SE)^ + {ts)^), where 
5E is the energy level change with Em = Eg + 6E 
and -Bill = Eg — 6E, Eg the ground state energy and 



SE cx (/^/yS - 0.5) Qualitatively, then, corre- 

sponds to the peak width of the envelope of Q-factor in 
Fig. El^d). Consequently, the stronger tg, the wider the 
range for the GHZ state. 

But the single-qubit tunnelling makes the coupled 
qubits unentangled. Thus we need to suppress the single- 
qubit tunnelling, while enhancing the three-qubit tun- 
nelling. In order to do so we need strong coupling as 
shown in Fig. 1 (d) , where the single-qubit tunnellings be- 
tween the ground and excited levels become suppressed. 
On the other hand the relevant parameter for tg is the 
Josephson couphng energy Eji. As Eji decreases, the 
barrier in the double-well potential in Fig. [T](d) becomes 
lower. It implies that, to get a larger value of tg, Eji 
should be smaller. But too small Eji makes some excited 
states, \saSbSc), unstable. We show the minimum Eji's 
in Table I for three representative E'/s. In fact, we found 
that for strong coupling case the excited states are stable 
for smaller Eji. For small E'j/Ej = 0.05 we obtained 
tf/Ej = 1.3 X IQ-^ and thus ig/i? « 5.4 x lO'^. But 
for larger value of E'j/Ej = 0.6 we obtained tg/t" ~ 6.4 
with ti/Ej = 8.0 X 10~^. Hence for strong coupling case 
we can expect higher Q-factor for the GHZ state. 

In Table I the peak widths for both GHZ and W states 
are calculated at 95% of the maximum value of Q-factor, 
which are approximately proportional to tg and t2, re- 
spectively. During the Rabi oscillations the fluctuation 
of flux is estimated to be in the order of 10~^[^o/ Hz^^^] 
[2^ and 1// critical current fluctuations of the Joseph- 
son junctions is rather weak. In recent experiments for 
flux qubits, the flux amplitudes are controlled up to the 
accuracy of 10^^<I>o. In this respect, the peak width, 
(5/^ - 4 X 10-^, for E'j = O.GEj wiU be sufficient to 
observe the GHZ state experimentally. 

W state.— In Fig. [11(b), we present the blue 
(light gray) triangle whose corners correspond to the 
three states consisting of the W state, |^'w) — 
(ITii) + liti) + liiT)) /V3. The blue (light gray) tri- 
angle intersects (j)~f-axis at 0^ > 0. Actually, there is 
another intersection plane with < for another possi- 
ble W state. For simplicity, we will focus on the W state 
on the blue (light gray) triangle plane. The effective po- 
tential at the plane of the blue (light gray) triangle for 
the three states is drawn in Fig. [Ije). Energetically, in 



E'j 


Eji 


,a 
H 


peak width 
GHZ {Sf-,) 




peak width 
W {5fa) 


0.05 
0.1 
0.6 


0.7 
0.75 
0.58 


7.0 xlO"'' 
2.6 xlO"'' 

5.1 xlO"^ 


~5 x 10"^ 
~2 X 10"** 
~4 X 10"* 


6.3 xlO"* 
1.0 xlO"* 



-10-'' 
~2 X 10"^ 




TABLE I: Peak widths for Q-factors of GHZ state in Fig. [g^d) 
and of W state in Fig. EJd) at 95% of the maximum values. 
Here the unit of E'j, Eji, and t is Ej. 
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our model, the energies of three states are higher than 
those of the two states consisting of the GHZ state, i.e., 
the ground state. Then, the W state can be observed in 
an excited state. 

Let us discuss how a W state can be reahzed in an 
excited state. At the coresonance point {.fa, fp, f-i) — 
(0, 0, •\/3/2), the six states except for , |tTT)} are 

degenerated in the second excited state. The six states 
are classified into two classes, {|TTi)7|TiT)iUTT)} with 
= 1/2 and {|ni),|in),|iiT)} with = -1/2. 
Hence, the two classes of the six states can be sepa- 
rated by applying an additional flux A/^. The three 
states of each class can form a W state. As shown in 
Fig. m^e), the two-qubit tunneling amplitude cre- 
ates the W-state, while the single-qubit tunneling t\ 
destroys the W-state because it induces a superposi- 
tion of states of the two classes. If other small tun- 
nelings are negligible, then, the Q-factor is given by 
g(|*)) ~ 8(1 -I- 2.5(tVA£^)^)/9(l + iit\/AE)^). From 
QQ'i'w)) = 8/9, it turns out that \AE\ oc \ fj/V3- 1/2| 
should be much larger than t\. Therefore, a sufficient 
A/^ are needed to generate a W state. Actually, we 
found that \Afj\ > 0.01 is sufficient to show the gen- 
eration of a W state (Figs. [S^a) and (b)). In Fig. 
|3Kb), the W state is formed slightly away from the point, 
/a = //3 = 0. For a relatively weak coupling the single- 
qubit tunneling as well as becomes larger. Thus, 
an additional small flux « 0.0004 will break the sym- 
metry so that the state j^f) closer to the W state would 
be formed. In Fig. ^c) we can see the W state around 

U = o. 

On the contrary to the GHZ state where the range 
of f.y is critical for experimental observation, for W state 
the range in (/□,, /^)-plane is important as shown in Figs. 
|3Ka) and (b). Fig. [3Kd) shows the Q-factor for W-states 
whose peak widths depend on the value of the two-qubit 
tunneling amplitude <2 (Table I). As Ej decreases, 
becomes larger. However, if the coupling strength be- 
comes too weak, the two classes with Sz = ±1/2 will 
become overlapped with each other through the single 
qubit tunnelling t\ so that the W state may readily be 
broken. Hence, as a consequence of compromise, the 
W state emerges for an intermediate coupling strength, 
Ej = Q.05Ej, with rather broader peak width as shown 
in Table I. 

Discussions and summary.— The quantification of en- 
tanglement can be done by using the state tomogra- 
phy measurement 0, [2^ . Recently for capacitively cou- 
pled phase qubits the tomography measurement has been 
done [91] , where they simultaneously measure the state of 
coupled qubits. For present coupling we expect that the 
similar tomography measurement can also be performed. 

The tripartite entanglement with superconducting 
qubits has not yet been achieved so far. For the bi- 
partite entanglement capacitively coupled phase qubits 
showed high fidelity in a recent experiment [9[, while for 




FIG. 3: (Color online.) The Q-factors of the second ex- 
cited state with f-, = ^/3/2 - 0.01 for (a) E'j/E.j = 0.1 and 
E.ji/E.j = 0.75 and (b) E'j/E.j = 0.05 and E.ji/E.j = 0.7. 
Note the different scales and positions of maximum Q-factor 
in both figures. Here the yellow (light gray) regions denote 
high Q-factors. (c) The coefficients of the eigenstate used 
in (b) are plotted along the dotted line with fp — —0.0004, 
which shows that W state is formed around fa = 0. (d) Cut 
view of Q-factors along dotted lines in (a) and (b) . The green 
dotted line indicates Q = 8/9 for W state. The peak width 
for Ej/E.j = 0.05 is much wider than that for Ej/E.j = 0.1. 



charge qubits only partial entanglement was observed. 
The interaction between phase qubits are XY-type inter- 
action which describes simultaneous two-qubit flipping 
processes. The two- or multiple-qubit tunnelling pro- 
cesses are essential for entanglement of qubits [l2j . How- 
ever, for charge qubits, the interaction are mainly Ising- 
type. We believe that this is the reason for weak en- 
tanglement in experiments with charge qubits. For three 
coupled phase qubits the lowest energy state is |000) state 
while the highest is |111) state. Hence the superposition 
between these two states will be negligibly weak, thus 
the GHZ state cannot be formed. But, since the other 
states, for example |100), |010), |001) states, are energet- 
ically degenerated, the W state could be obtained. 

In summary, we investigate a three superconducting 
flux qubit system. The GHZ and W states can be re- 
alizable in the eigenstates of the macroscopic quantum 
system. We show that while the GHZ state needs strong 
coupling strength, the W state can be formed at an opti- 
mized coupling strength. Moreover, to keep the tripartite 
entangled states robust against external flux fluctuations 
for feasible experimental realizations, the three coupled 
qubit system can provide relatively large three-qubit and 
two-qubit tunneling amplitudes for GHZ and W states, 
respectively. 

SYC acknowledges the support from the Australian 
Research Council. 
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